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ABSTRACT 

The Fermi LAT discovery that classical novae produce > 100 MeV gamma-rays estab¬ 
lishes that shocks and relativistic particle acceleration are key features of these events. These 
shocks are likely to be radiative due to the high densities of the nova ejecta at early times 
coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are ab¬ 
sorbed by neutral gas and reprocessed into optical emission, similar to Type Iln (interacting) 
supernovae. Gamma-rays are produced by collisions between relativistic protons with the 
nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic 
electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic 
particle energy into LAT gamma-rays is at most a few tens of per cent. The measured ratio 
of gamma-ray and optical luminosities, Ly/Lopt, thus sets a lower limit on the fraction of the 
shock power used to accelerate relativistic particles, fnth- The measured value of Ly/Lopt for 
two classical novae, V1324 Sco and V339 Del, constrains e„ih > 10“^ and > 10“^, respec¬ 
tively. Leptonic models for the gamma-ray emission are disfavored given the low electron 
acceleration efficiency, enth ~ 10 - 10“^, inferred from observations of Galactic cosmic 

rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(enth/0.01) ' and 
> 10(enth/0.01)“' per cent of the optical luminosity is powered by shocks in nova Sco and 
nova Del, respectively. Such high fractions challenge standard models that instead attribute 
all nova optical emission to the direct outwards transport of thermal energy released near the 
white dwarf surface. We predict hard ~ 10-100 keV X-ray emission coincident with the LAT 
emission, which should be detectable by NuSTAR or ASTRO-H, even at times when softer 
< 10 keV emission is absorbed by neutral gas ahead of the shocks. 
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1 INTRODUCTION 


As “guest stars” to our ancestors, novae have been observed since 
antiquity (|Zwicky||1936|. Powered by runaway nuclear burning 
on the surface of white dwarfs (e.g. [Starrheld et al.|[T972| l, nova 
eruptions are usually accompanied by the expulsion of matter out¬ 
wards at high velocities >10^ km s“* ( Shore|2013 and references 
therein). Most of the radiated energy in a nova outburst occurs at 
optical and ultraviolet wavelengths and is generally attributed to 
the outwards transport of thermal energy produced near the white 
dwarf surface ( |PriWik|[l986[ [Prialnik & Kovetz||1995[ |Starrfield] 
|et al.|20()0t . Dynamical stellar evolution and hydrodynamic calcu¬ 
lations are used to interpret basic features of nova light curves, such 
as how their rate of evolution depends on the white dwarf mass. 


* E-mail: bmetzger@phys.columbia.edu 


central temperature, and accretion rate (e.g., |Yaron et al.||2005[ 
|Hillman et al.|[20T4l l. Despite notable successes, models of nova 
outbursts remain plagued by uncertainties such as the efficiency of 
convective mixing and the assumed prescription for mass loss. This 
makes it challenging to accurately predict the quantity, velocity, and 
time history of matter ejected from the white dwarf surface. 

The Large Area Telescope (LAT) on the Fermi satellite re¬ 
cently discovered that classical novae emit > 100 MeV gamma- 
rays, occurring at times coincident to within a few days of the op¬ 
tical peak and lasting for several weeks ( [Ackermann et al.|2014| ). 
Luminous continuum gamma-ray emission indicates the presence 
of strong shocks which accelerate particles to relativistic energies, 
producing gamma-rays via leptonic or hadronic processes (e.g., 
Tatischeff & Hernanz||2007^ . Radio VLBI imaging by |Chomiuk] 
et al.| ( |2014| > recently detected compact, high brightness tempera¬ 
ture knots from V959 Mon, confirming the presence of non-thermal 
shock emission in a gamma-ray producing nova. 
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Evidence for shocks in novae has existed well before Fermi 
in the form of hard X-ray emission (e.g., |Mukai & Ishida||2001| 
|Mukai et al.||2008| l, early peaks in the radio emission which are 
inconsistent with thermal emission from freely expanding photo- 
ionized gas (e.g., |Taylor et al.||1987| |Krauss et al.||2011| |Weston| 
|et al.|2013t , and coronal line emission during the nebular phase re¬ 
quiring a harder source of ionizing radiation than expected from 
the cooling white dwarf (e.g., [Shields & Ferland|197^ . However, 
since most of these signatures are observed on timescales of months 
or later after the eruption, shocks appear to have been relegated to 
a mere side feature of the main thermonuclear event. The Fermi 
discovery of luminous shocks in coincidence with the optical peak 
unambiguously establishes their importance to the qualitative pic¬ 
ture of nova eruptions. The dearth of X-ray and radio shock sig¬ 
natures at times coincident with the gamma-rays is unsurprising in 
retrospect due to the high bound-free and free-free optical depths 
created, respectively, by large columns of neutral and ionized gas 
at early times jHillman et al.|2014{ [Metzger et al.|20r4) . 

Shocks were unexpected in classical novae because the pre¬ 
eruption environment surrounding the white dwarf is occupied only 
by the low density wind of the main sequence companion star, re¬ 
quiring a different source of matter into which the nova outflow 
collides. One physical picture, consistent with both optical (e.g., 
[Schaefer et al.[201^ and radio imaging (e.g. |Chomiuk et al.[2014[(, 
and the evolution of optical spectral lines (e.g., [Ribeiro et al.|2013[ 
[Shore et al.[[2013) , is that the thermonuclear runaway is first ac¬ 
companied by a slow ejection of mass with a toroidal geometry, 
the shape of which may be influenced by the binary companion 
(e.g.,[Livio et al.[1990|[Lloyd et al.[1997|l. This slow outflow is then 
followed by a second ejection or a continuous wind (e.g., [Bath &[ 
[Shaviv[1976^ with a higher velocity and more spherical geometry. 

Assuming the expanding ejecta cools sufficiently, the subse¬ 
quent collision between the fast and slow components produces 
strong “internal” shocks within the ejecta which are concentrated in 
the equatorial plane. The fast component continues to expand freely 
along the polar direction, creating a bipolar morphology (Fig.[TJ. 
Interestingly, modeling of the symbiotic gamma-ray novae V407 
Cyg ( [Abdo et al.|2010[ ) also appeared to require the presence of a 
dense equatorial torus ( [Martin & Dubus|2013), si milar to that in¬ 
ferred in classical novae (see also [Sokoloski et al.|20()8l [Orlando] 
[et al.[[200^ [Drake et al.[|2009l [Orlando & Drake1[2012[ for other 
evidence for bipolar ejecta in symbiotic novae). 

[Metzger et al.[ ( [2014^ present a semi-analytic model for nova 
shocks and their thermal emission, which they fit to the radio and 
X-ray data of the gamma-ray nova V1324 Sco. A key finding was 
that the shocks responsible for the radio maximum seen months af¬ 
ter the eruption could, at smaller radii and earlier times, also power 
the nova optical emission. One can generalize the basic argument 
as follows. Gas ahead of a shock with power sufficient to create the 
observed gamma-ray emission is necessarily dense. Most of the ki¬ 
netic energy dissipated by shocks moving through a dense medium 
is radiated as thermal X-rays (a so-called “radiative” shock). These 
X-rays are absorbed by neutral gas ahead or behind the shock, re¬ 
processing their energies to lower, optical frequencies, where the 
lower opacity readily allows their energy to escape. The observed 
gamma-ray luminosity is typically only a fraction ~ 10“"^ - 10“^ 
of that emitted as optical/UV radiation (^. However, since only a 
fraction of the shock power is used to accelerate relativistic parti¬ 
cles and only a fraction of that is radiated in the FAT bandpass, one 
concludes that a significant fraction of the nova optical emission is 
shock-powered. 

The above result also implies that combined optical and 


gamma-ray data from novae can be used to probe particle accel¬ 
eration at non-relativistic shocks. The basic concepts of diffusive 
shock acceleration (DSA) were developed almost forty years ago 
(e.g., [Axford et al.p^ [Bell|1978[ [Blandford & Ostriker|1978) . 

However, only recently have plasma kinetic simulations seen the 
self-consistent development of the DSA cycle in collision-less non- 
relativistic shocks (e.g.,[Caprioli & Spitkovsky|2014a[[Kato|2014[ 
[Park et al.|2^ . As we shall discuss, gamma-ray novae provide 
a real-time laboratory for testing the predictions of DSA theory in 
ways complementary to traditional methods, such as the modeling 
of supernova remnant emission. 

This paper solidifies the above theoretical arguments and puts 
them into practice using data from the classical novae VI324 Sco 
and V339 Del (hereafter, ‘nova Sco’ and ‘nova Del’, respectively). 
Though well-studied at many wavelengths, we do not consider the 
nova V959 Mon because of the lack of early optical coverage con¬ 
current with the gamma-ray detections. The salient properties of 
nova shocks are reviewed in ^ including their radiative nature 
the high efficiency with which shock power is radiated at 
optical/UV frequencies ( ^2.2[ (, and the efficiency of particle accel¬ 
eration and gamma-ray production ( |2.3^ . A combined analysis of 
the optical and gamma-ray data of VI324 Sco and V3229 Del is 
presented in ^ using our theoretical framework to constrain the 
minimum fraction of the shock power placed into relativistic parti¬ 
cles and the minimum fraction of nova optical light curves that are 
shock powered. In ^we discuss our results and their implications 
for particle acceleration at non-relativistic shocks and for the power 
source behind nova optical light curves. 


2 SHOCKS IN GAMMA-RAY NOVAE 
2.1 The shocks are probably radiative 

The slow outflow from a nova with velocity V4 = lO^Vg km s“' 
expands to a radius 


Rq = V4? ~ 6 X lO'^twkVg cm 


( 1 ) 


by a time t = week. The density in the slow ejecta of assumed 
thickness ~ Rej and hydrogen-dominated composition is given by 


Mg 


~9x cm-\ 


( 2 ) 


where /ao ~ 0.5 is the fraction of the total solid-angle subtended 
by the outflow (Fig.j^ and is the ejecta mass, 

normalized to a value characteristic of those measured by late ther¬ 
mal radio emission ( [Hjellming et al.[ri979[ [Seaquist et al.[[1980^ . 
We assume a radially-thick outflow because the resulting low den¬ 
sity makes this the most conservative assumption when arguing for 
radiative shocks. A wide solid angle /ao ~ 0.5 is also the most 
conservative assumption when arguing for radiative shocks and is 
needed to produce a high efficiency for converting the kinetic en¬ 
ergy of the outflow to shock, and hence gamma-ray, luminosity. 

A faster outflow (hereafter “wind”) of mass loss rate M and 
velocity Vi ~ 2v4 collides with the ejecta from behind. The relative 
velocity between the fast and slow outflows is unlikely to greatly 
exceed this assumed value ~ V4 ~ 10^ km s“', but if it is much 
lower than assumed this would only strengthen our subsequent ar¬ 
gument that the shocks are radiative, since the cooling time is an 
increasing function of the shock velocity. The density of the wind 
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Figure 1. Schematic diagram of the proposed geometry (side view) of shock 
interaction in classical nova outflows (see text). A slow outflow with veloc¬ 
ity V4 < 10^ km s“^ is ejected first, its geometry shaped into an equatorially- 
concentrated torus, possibly due to interaction with the binary companion 
of the white dwarf (e.g., |Livio et al.|199Q^ . This outflow is followed within 
a few days by a faster outflow or continuous wind with a higher velocity 
V] ~ 2 v 4. The fast and slow components collide, produce a forward-reverse 
shock structure. Kinetic energy dissipated by the shocks is radiated as ther¬ 
mal X-rays, which are absorbed by neutral gas ahead or behind the shock 
and re-radiated as thermal optical/UV emission ^Metzger et al.|2014^ . A 
small fraction, e^th 1, of the shock power is used to accelerate non- 
thermal ions or electrons, which radiate gamma-rays by interacting with 
ambient gas or radiation, respectively. 


at the collision radius (~ radius of the slow ejecta) is given by 

~ , Jf - 2 X (3) 

AnRljripVi 

where M = 10 “^M_ 5 Mq wk“' is normalized to a value resulting in 
the ejection of ~ lO^^M© over a week. For instance, a total mass 
4 X IO^^Mq was ejected in the “fast” component of V959 Mon 
( [Chomiuk et al.|[2014| l. Our assumption that the fast outflow is a 
steady wind is also the most conservative one when arguing for 
radiative shocks; for a fixed fast ejecta mass, the wind density at 
the reverse shock is smaller than if the ejection had occurred over a 
shorter duration. 

The interaction drives a forward shock (FS) through the slow 
shell and a reverse shock (RS) back through the wind (Fig.[T](. As¬ 
suming the shocks are radiative, the post shock material is com¬ 
pressed and piles up in a central cold shell sandwiched by (the 
ram pressure of) the two shocks. The FS propagates at a veloc¬ 
ity Vf = Vc - V4 <K V4, while the velocity of the reverse shock is 
Vc is the velocity of the cold central shell (see 
although be wary of notationaQ differences). 


Vr = Vi - Vc, where 
Metzger et al.|2014 


* In particulai, |Metzger et al.|j2014f define Vsh as the velocity of the central 
shell (Vc in this paper). In this paper Vsi, instead denotes the velocity of the 
shocks. 


Flere we have assumed the shocks are radiative, for which the ve¬ 
locity of the immediate post shock gas slows to that of the cold cen¬ 
tral shell, such that the shock velocity equals the difference between 
the velocity of the upstream gas and the central shell. Hereafter, the 
velocities of both shocks are parametrized as Vsh = ? 7 V 4 , where we 
expect that rj < \ for the FS and 77 ~ 1 for RS. Also note that for pa¬ 
rameters of interest, the amount of kinetic energy dissipated by the 
RS is greater than that by the FS due to the higher velocity of the 
RS, potentially favoring the RS as the site of particle acceleration. 

The shocks heat the gas to a temperature 

Tsh - ~ 1.4 X lO’v^T,^ K (4) 

iD/C 

and compresses it to a density = 4nej (FS) or 4n„ (RS), where 
we have taken = 0.62. 

Gas cools behind the shock on a characteristic timescale 
_ 3kT,J2 _ ( 1.6x10V^M:]4s FS, 
„,,A(r,h)~l 7.8 x10XM_14s RS, 


where A = AoT*^^ = 2 x erg cm^ s“' is the cooling 

function and we have assumed 77 = 1 in the RS case. Our cool¬ 
ing function includes just free-free emission, which is conserva¬ 
tive from the standpoint of radiative shocks because line cooling 
contributes a comparable or greater cooling rate at temperatures 
< 3 X 10 ’ K (e.g. Schure et al.|2009V Here we have assumed that 


electrons and protons are efhciently coupled behind the shock, as 
is justified because the timescale for Coulo mb energy e xchange 


lAXlf Irish s (NRL Plasma Formulary; Huba|20o'^ is short 


compared to the expansion timescale lexp ~ Rej/v 4 ~ t. 


3 X 10-'M:>«7734 FS, 
t ~\ 1.2x RS, 


( 6 ) 


where we have (conservatively) assumed a pure hydrogen compo¬ 
sition. 

Whether the shocks are radiative depends on the ratio of the 
cooling timescale to the expansion timescale. 


f 2.7xl0-VX-i4 FS, 

^ f ~ \ 0.13v^M:^f„k RS, ^ ’ 

For the characteristic range of parameters M _4 ~ 0.3 - 3 ([Seaquist] 
& Bode|2008|, vg < 3, 77 < 1 and timescales corresponding to the 
gamma-ray emission t < few weeks, one concludes that the FS is 
likely to be radiative ix < 1). It is less clear that the RS is radiative, 
since in principle for vg > 2 and a very low wind mass loss rate 
M _5 < 1 one can have;y > 1. 

The argument for radiative shocks can be strengthened by con¬ 
sidering additional constraints. First, the optical depth of the outer 
unshocked ejecta at optical frequencies, 

fopt ~ ^yj^ejRej^opt ~ 0.9A/_4Vg (8) 

must exceed unity, where /Cop, ~ 0.1 cm^ g“* is the optical opacity 
set by Doppler-broadened iron lines ^Pinto & Eastman|2000| ). This 
constraint is motivated by the lack of clear evidence, e.g. emission 
lines, for hot, energetic shocks in optically thin regions during the 
optical peaks of novae and during most epochs of gamma-ray de¬ 
tection. Over the first few weeks or longer following outburst (i.e., 
prior to the drop of the “iron curtain”), nova spectra are character¬ 
ized by broad P-Cygni line profiles, indicating the presence of a 
pseudo-photosphere. 

Equation l|7j can be rewritten in terms of T(,pt as 

9'"?/"^F277v^„7ropt ^ r 2 X 10-377v2T-pV FS, 
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For the forward shock to occur helow the photosphere (Topt > 1) 
yet not be radiative Cy > 1 ) would thus require unphysically high 
ejecta velocities, V 4 > 10,000 km s“'. The RS is also radiative if 
T(,pt > 1 on timescales of a couple weeks for M _4 < 3 as long as 
M_5 > 0.2. 

A second requirement is that power dissipated by the shocks 
Lsh = 97rR^jUsh'npV^p/32 exceed the observed gamma-ray luminos¬ 
ity Ly by an efficiency factor l/enthty > 100 accounting for the 
fraction of the shock power radiated as gamma-rays, where ac¬ 
counts for the efficiency of the shock in accelerating non-thermal 
particles and Sy accounts for the energy of the latter radiated in the 
LAT bandpass (see & Equation 0 can be written in terms of 
Tsh “ Ty/6nthfy ^tS 


A 


81;r ej ^nth 


256 V3 

0.0277S'gfwk 


LyAo 

Ly I 6,1th j 

1036 erg s-‘ ) \ 0.01 /' 


( 10 ) 


For measured values Ly ~ 3x lO^^ erg s“* i Ackermann et al. 2014| 
the shocks are thus radiative Oc < 1 ) on a timescale of weeks for 
6 nth 6 y < 0.01 if the velocity of the shocks is < 2,000 km s“'. 

In conclusion, both forward and reverse shocks are likely to be 
radiative at times corresponding to the observed gamma-ray emis¬ 
sion, although this statement is the most secure at the earliest times. 
One can nevertheless keep the following arguments fully general by 
introducing the shock radiative efficiency/rad = (l-l-5^/2)“*, which 
equals unity for x ^ 1 but scales oc 1/x for x s> 1 (|Metzger et al.| 
|2^ . 


2.2 The shock power will mostly emerge as optical radiation 

Absent the immediate presence of a shock, the bull0of the nova 
ejecta is neutral at early times because the timescale for radia¬ 
tive recombination, tree ~ l/"cjttrec ~ ^ 1 ® 

tremely short compared to the evolution timescale ~ weeks, where 
ttrec ~ 10^*323 cm3 ^-1 approximate radiative recombination 
rate for hydrogen-like species of charge Z ( [Osterbrock & Ferland| 
2006 1 . X-rays of temperature < 10’ K (eq. |4|) are thus ab¬ 
sorbed by high columns of neutral gas ahead or bMiind the shock, 
before being re-radiated as line emission at lower frequencies. If 
the shock is radiative then the timescale for reprocessing to lower 
temperatures is also necessarily short because the line cooling func¬ 
tion A(T) increases rapidly with decreasing temperature down to 
T ~ 10“* K (e.g. Schure et al.|2009| . However, harder X-rays with 
energies > 10 keV should be free to escape even at these early 
times due the decreasing bound-free cross section at high photon 
frequencies (^. 

Radiation escapes efficiently once it reaches the optical/near- 
UV band due to the much lower opacity at these frequencies com¬ 
pared to the UV/X-ray opacity of neutral gas. Reprocessed optical 
radiation will furthermore not be degraded by PdV losses provided 
that the photon diffusion timescale ~ Rej^opt/c is shorter than 
the expansion timescale texp ~ f^ej/t '4 ~ L where Top, is the opti¬ 
cal depth (eq. j^). Equating these two, optical radiation escapes 


without adiabatic losses after a time jAmett|1982|> 


top, « 


( 11 ) 


The timescale top, also sets the minimum timescale for the optical 
lightcurve to rise after the onset of the explosion. 

The onset of gamma-ray emission from V1324 Sco and V339 
Del was delayed with respect to the peak of the optical radiation 
by a few days i Ackermann et al.|2014 1 . The dominant proces^by 
which gamma-rays of energy 6 y ~ 0.1 - 3 GeV are attenuated is 
inelastic electron scattering, for which the Klein-Nishina opacity 
is /Tie ~ 10“3(f^/GeV)“*k-es ~ 4 X 10“3(f^/GeV)“*k-op,. The ejecta 
will thus remain opaque to gamma-rays of energy gy until Top, < 
250(ey/GeV), as occurs after a time ty which exceeds that of the 
nominal optical rise time (eq. |1H) by the ratio 


Gpt \ t'w^op. 


1/2 


■■(civ) 


- 1/2 


-1/2 


( 12 ) 


i.e. ty ~ couple days for ~0.1 GeV. 

From the above we can draw two key conclusions: (i) a shock 
that produces gamma-ray emission which is not absorbed and 
hence observable (gamma-ray optical depth Ty < 1) necessarily 
radiates the bulk of its dissipated thermal energy without adiabatic 
losses at optical frequencies (Top, < C/V 4 ). In other words, at times 
when the column of ejecta is sufficiently low to allow gamma-rays 
to escape unattenuated, it is also sufficiently low to allow ther¬ 
mal radiation to escape without significant losses to PdV work, 
(ii) gamma-ray absorption (t^ > 1 ) at early times may explain the 
delayed onset of the gamma-ray emission. If true, the latter im¬ 
plies that the optical emission near peak can be shock-powered, 
even if gamma-ray emission is suppressed at this time. Alterna¬ 
tively, the gamma-ray delay may result from the finite timescale re¬ 
quired to accelerate the gamma-ray producing particles at the shock 
(G. Dubus, private communication). 


2.3 Partitioning the shock energy 


A large fraction /ad ~ 1 of the total power Ljh dissipated by shocks 
goes into thermal X-rays, which are absorbed and re-radiated as 
optical emission ( §2.2| l. A much smaller fraction 6 „,h ^ 1 goes into 
accelerating non-thermal ions or electrons. The fraction of non- 
thermal power 6 nth radiated as gamma-rays, Ly = 6 n,heyLsh, depends 
also on a factor 6 ^ < 1 accounting for the radiative efficiency of the 
accelerated particles and the fraction of the total gamma-ray emis¬ 
sion emitted in the LAT bandpass. Combining these expressions, 
the fraction of the total nova optical luminosity powered by shocks 
can be written 

_ Lopt.sh /rad Ly 

M=—. - ~ - f —. (13) 

^opt ^y^nth ^opt 


Once ey is specified based on the assumed emission process 
(hadronic or leptonic), the observed ratio LylL^^x sets a lower limit 
on the value of enth» i-®- 


1 Ly 

fs\i ^ ^nth ^ ^nth.min ~ ? 

^y ^opt 


(14) 


assuming a radiative shock (/rad = 1). Alternatively, if the value of 
6 nth is assumed, then the measured value of Ly/Lopt determines /sh. 


^ Although the bulk of the ejecta is neutral at early times. X-rays from the 
shocks fully ionize a thin layer of gas just ahead of the shocks jMetzgeret al.| 
|2Q14) . Thus, we do not expect the shock to be modified by the presence of a 
neutral upstream medium, as may occur in some supernova remnants (e.g., 
|Blasi et al.|2Q12| . 


^ Photo-pion production is the dominant source of opacity for higher en¬ 
ergy gamma-rays 6y > 3 GeV, for which the opacity is ^ 3 x 10“‘^/('es 
(e.g. [Anchordogui et al.|20Q2||Montanet et al.|1994) . Photon-photon pair 
creation opacity becomes important at TeV energies (eq. 1251). 
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Gamma-ray novae as probes of particle acceleration 5 


We now consider what values of and 6 y are expected in hadronic 
and leptonic scenarios, respectively. 


2.3.1 Hadronic scenario 


The momentum distribution f{p) oc of non-thermal protons ac¬ 
celerated via Diffusive Shock Acceleration (e.g., |Blandford & Os-| 
|triker|1978| l corresponds to an energy distribution 

‘^El oc I ~ . (15) 

dEp ' ( constant, iripC «Ep < iimax. 


that concentrates most of the non-thermal energy in relativistic par¬ 
ticles. The |Hill^jl984| > criterion sets an upper limit on the maxi¬ 
mum proton energy 

- VX ,16) 

where Ssh is the magnetic field strength near the shock, which is 
estimated by assuming that magnetic pressure fis|,/ 8 ;r is a fraction 
fB of the thermal pressure of the post-shock gas. 

A reasonably high value of eg is expected if ion acceleration is 
efficient, because cosmic-ray induced instabilities amplify the mag¬ 
netic field for some distance behind the shock to a significant frac¬ 
tion of its equipartition value. [Caprioli & Spitkovsl^ ( |2014b] l show 
that the magnetic field near the shock is amplified from its initial 
upstream value So, according to (Ssh/So)^ ~ SCnthMA (their eq. 2 ), 
where = Vsh/v^ is the Alfvenic Mach number of the upstreanj^ 
flow in the shock frame, Va = So/ Apnrnfnl, and is the fraction 
of the shock dissipated kinetic energy placed into relativistic ions. 
The equipartition fraction of the shock region can thus be written 
as 

’ ( 17 ) 

where eg^o = E^KAnnumpV^) is the ratio of magnetic energy to ki¬ 
netic energy within the initial, unshocked outflow. 

The value of eg o within the unshocked nova ejecta can be con¬ 
strained as follows. Assume that the magnetic energy density of the 
outflow Ub is a fraction eg,* of the bulk kinetic energy oc pv ^/2 
near the stellar surface (r = R^,), where p and v are the outflow den¬ 
sity and velocity, respectively. Flux-freezing within the 3D expand¬ 
ing flow (p oc r“3) dilutes the field strength according to S oc 1/r^, 
such that UbIU]^ = ee.o at the radius of the shock ~ will be at 
most a factor of f?*//?ej smaller than 6g,*. For characteristic values 
of R^j ~ lO*'^ cm and the initial radius R^, > 10 ^ cm of the out¬ 
flow set by the white dwarf surface, we expect that eg,o > 10 “^eg,*. 
Thus, for Cnth = 0.1 and rj = 1, even an extremely small value of 
eg * > 10 “*® appears sufficient to accelerate protons to a maximum 
energy £,^ 3 * > 10*** eV (eq. |16|) large enough to explain the high¬ 
est LAT-detected photon energies via no creation (|Ackermann et al.| 
|2m4l l. 

We caution, however, that the bulk of the ejecta is neutral and 
hence may not be able to support the strong turbulent magnetic field 
needed to accelerate particles. In this case the radius R^j entering the 
Hillas criterion in equation should be replaced by the width of 
the X-ray ionized layer ahead of the shock. [Metzger et al.| j2014^ 
estimate the latter to be A‘™ ~ 9xlO*’;7V8(nej/10*** cm“-’)“* cm based 


'* For concreteness, we focus here on amplification of the magnetic field of 
the slow ejecta at the forward shock, although similar considerations apply 
to amplification of the field in the fast outflow at the reverse shock. 


on ionization by the free-free emission (their eq. [45]). Replacing 
with Ajon in equation l|16|l results in 


^Bsht^shAic 


1.0 X 10‘*eV;7® 


(V 

( '1 

\10-6/ 

\ 10*® cm“® / 




Thus, for typical values of - lO’ - 10*** cm ® (eq. j^, j^) 
and vg ~ 1 , one finds that higher field strengths of eg > 10 ““ - 
10“"* (eg,* > 10“^ - 1) is now required to have £max > 10*** eV. 
Note, however, that jMetzger et al.]j2014^ do not take into account 
additional ionization of hydrogen by the UV radiation produced by 
line cooling or by the reprocessed X-rays, which could be important 
at early times of relevance because the ejecta is opaque to X-rays. 
This could increase the energy of the ionizing radiation, and hence 
the thickness of the ionized layer and £miix, by up to a factor of 
~ kTshl(2 Ryd) ~ 5t)Tfv\, where Ryd = 13.6 eV is the characteristic 
ionization threshold energy for hydrogen, relaxing constraints on 


The fraction of the shock power placed into non-thermal ions 
fnth can be estimated from observations of other non-relativistic 
shocks, such as gamma-ray emission in supernova remnants 
(e.g. [Ack ermann et al.||20l4l l. In Tycho, for instance, [Morlino 
[Caprioli j20I2) infer e^h ~ 0.1 if the gamma-rays are hadronic 
in origin. Hybrid (kinetic ions - fluid electron) simulations of non- 
relativistic shocks also find fnth ~ 0.1 - 0.2 for cases in which the up¬ 
stream magnetic field is quasi-parallel to the shock normal ([Capri-| 
[oli & Spitkovsky|2014a) . However, if the nova ejecta is character¬ 
ized by a phase of approximately steady-state outflow, then radial 
expansion in two dimensions will produce a magnetic field domi¬ 
nated by its toroidal component, i.e. perpendicular to the outflow 
and hence shock direction (Fig. [^. For such quasi-perpendicular 
magnetic field geometries, [Caprioli & Spitkovsky| ( [2014a[ l infer 
much lower proton acceleration efficiencies (consistent with zero), 
a point we return to in 0 

Relativistic protons accelerated at the shocks produce pions by 
colliding with effectively stationary protons in the ejecta. The mean 
time between interactions is giv en by tp_p = {Kp^pnip n^py^, where 
Kp-p ~ 0.025 cm“^ g“* X! Kopt!^ I Kamae et al.|2006 1 is the opacity 
for inelastic proton collisions of energy Ep > GeV. The number of 
collisions a proton experiences over an expansion timescale f/tp_p 
thus exceeds unity until after a time 


?p-p ~ 8.2 weeks. (19) 


At times t « tp_p, protons lose their energy to pion production 
instead of adiabatic expansion; the ejecta thus acts as an efficient 
“calorimeter” for converting relativistic protons into gamma-rays. 

The fraction of p — p collisions producing gamma-rays is 
the product of the fraction of inelastic interactions, cr^^lfa^^ + cr^) ~ 
0.5 - 0.8 (where cr^ and Cie are the elastic and inelastic cross sec¬ 
tions, respectively; [Kamae et al.|2006) , and the fraction =: 1/3 of 
inelastic events placed into the p + p^n^^y + y channel. We 
thus typically expect that Cy ~ 0.2-0.3, depending on the details of 
the accelerated proton spectrum, while a lower value could result if 
protons are not efficiently trapped, i.e. at times t > tp_p. 

To summarize, we expect enthfy < 0.03 in hadronic scenar¬ 
ios if the magnetic field within the ejecta is perpendicular to the 
shock plane, but this value may be considerably lower if the field 
is instead parallel to the shock plane, as would be expected given a 
phase of quasi-steady outflow from the white dwarf surface. 
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2.3.2 Leptonic scenario 


In supernova shocks, the ratio of the energy placed into rela¬ 
tivistic protons to that in relativistic electrons is estimated to be 
Alep ~ 10^“* - 10“^ based on observations of individual remnants 
( |Vdlk et al.|2005[|Morlino & Caprioli|2012^ and synch rotron emis¬ 
sion from Galactic cosmic rays ( |Beck & Kjause|2005[[Strong et ^ 
|2010| l. This value is consistent with recent particle-in-cell simula¬ 
tions of non-relativistic shocks, which find /^ep ~ 10“^ when extrap¬ 
olated to shock velocities v^h/c < 0.01 characteristic of those in no¬ 
vae ([Park et al.|2014[[Kato|2014l>. Assuming shocks accelerate pro¬ 


tons with a maximum efficiency ~ 0.1 (5 2.3.1 


corre- 

3 


sponds to a relativistic non-thermal electron fraction of e„th <10 
To Compton upscatter optical seed photons of energy E^pt ~ 
eV energy to ~ 0.1 - 10 GeV requires electrons with energy E^, = 
and Lorentz factor y,, ~ 10“* -10^. The ratio of the Compton 
cooling timescale of the electron to the expansion timescale is given 
by 


tcool.IC 


2>mgC 


t ^a-TUoplJet 


^ 0 . 07 t- 


(7e r 

* / Lopt \ 

uo*/ 

\10^®ergs“*/ 


2 

Vgtwk; 


( 20 ) 

where 1/opt ~ Ft = LoptTopt/4;rcA^j is the radiation energy den¬ 
sity near the shock (for shocks below the photosphere, Topt > 1), 
where we have used the fact that the conserved radiative flux 
F oc dll apt I dr in the diffusion approximation implies that Ft is 
constant from the photosphere (t ~ 1) to the shock depth t = Topt. 
Equation i20i shows that y^ > 10“* electrons are in the fast cool¬ 
ing regime (rcooi hxp) for luminosities Lopt > 10“*^ erg s“' and 
timescales t ~ weeks of relevance. 

Relativistic Bremsstrahlung emission represents an alterna¬ 
tive, and possibly dominant, leptonic emission mechanism. The 
cooling rate of a relativistic electron of Lorentz factor y^ of f protons 


in the nova ejecta is given by = 3 x 10 erg s ' i Rybicki 

|& Lightman|1979] l, resulting in a cooling time 

^2 / „ \-l/2 


tcool,rb ygUlgC 


: 5.1 X 10 " 




( 21 ) 


t 4i-b ' 10“* 

where we have used equation A comparison between equa¬ 
tions ( |20^ and l |20| l shows that relativistic bremsstrahlung emission 
can exceed that from Inverse Compton scattering for typical values 
of the nova luminosity, outflow velocity, and ejecta mass. Regard¬ 
less of whether Inverse Compton or bremsstrahlung emission dom¬ 
inates, the shock energy placed into electrons with the necessary 
Lorentz factors to produce the observed gamma-ray emission will 
in fact be radiated with high efficiency. 

Compton scattering of electrons with an accelerated spectrum 
dNeldEe oc Eg’’ , where E^ = yem^c^, produces a gamma-ray spec¬ 
trum 


dNy 

~dE 


£-{p+l)/2 

p;-{p+2)/2 


£ « £c 

E » £„, 


( 22 ) 


where £c ~ £optTc 7c is the minimum electron Lorentz fac¬ 
tor obeying the fast-cooling condition (tcooi.ic < Lxp; eq. j20|). 
A similar cooling energy can be defined in the case of relativis¬ 
tic brehmsstrahlung emission, below which the photon spectrum 
dNyjdE oc E-(2p-i)l2^ j g jjjg same as for Inverse Compton in the 
p = 2 case of greatest relevance (see below). 

The fraction of the total energy in relativistic electrons radi¬ 
ated in the LAT bandpass is given by 


JlO-t 

f'°' •EkEdE 
Jl dEe 


108-4p _ IQIO-SP 0.2, p > 
0.2 p = 


3) 


where we have assumed that that particles with 10“* < 7e < 10^ 
are the only ones contributing into the LAT bandpass and that they 
radiative with 100 percent efficiency, i.e. tcooi.ic « t or fcootrb t. 
Power-law fits to the gamma-ray spectra of three classical novae 
yield best-fit photon indices T ~ 2 - 2.5 ( |Ackermann et al.|2014| l 
corresponding to p = 2(r - 1) ~ 2-3, for which we estimate 
Cy ~ 0.2 - 10^“*. To be conservative we hereafter assume p = 2, 
corresponding to ey = 0.2. 

In summary, adopting fnth ^ 10"^ - 10"^ motivated by simula¬ 
tions and observations, we estimate that e„tbG ~ 10"** - 10"“* in the 
leptonic scenario. 


3 DATA 

This section describes our analysis of the gamma-ray and optical 
light curves of V1324 Sco and V3229 Del, events chosen as cur¬ 
rently being the only classical novae with both published gamma- 
ray data and contemporaneous optical coverage. We do not con¬ 
sider symbiotic novae in this analysis because the argument for the 
shocks being radiative is less secure given that the latter propagate 
into the extended wind of the giant companion star, instead of the 
potentially denser slow ejecta in the case of internal shocks. Our 
goal is to use the measured ratio of gamma-ray and optical fluxes 
to constrain the particle acceleration efficiency e„th and the fraction 
of the classical nova optical light curve powered by shocks /sh, fol¬ 
lowing the arguments outlined in §2.3| 


3.1 Optical 


For each nova we use photometric measurements at both opti¬ 
cal and near-infrared (NIR) wavelengths. Optical photometry was 
taken from the database of the American Association of Variable 
Star Observers ( |Henden|2013^ . Both nova Sco and Del had BVR 
band measurements, but only Del included I band. To supplement 
the NIR for nova Sco, we used JHK measurements from the Small 
and Moderate Aperture Research Telescope (F. M. Walter, private 
communication). Only photometric measurements coincident with 
the published Fermi gamma-ray light curve were used; specifi¬ 
cally August 16-September 14 2013 and June 15-June 30 2012 
for V339 Del and V1324 Sco, respectively. Photometric reddening 
corrections were applied to the datasets using an E{B - V) = 0.2 
for V339 Del ( |Munari et al.l20T3l l and E(B - V) = 1.0 for V1324 
Sco (|Finzell et al.|2015f. Filter specific reddening corrections for 
the optical/NIR were taken from|Schlafly et al.|(|2011) assuming 
Rv = 3.1. 


The total optical/NIR flux each night, was determined by 
approximating the SED as a blackbody, fitting a temperature and 
normalization, and integrating over frequency. The early-time flux 
obtained for V339 Del are consistent with the values found in the 
more detailed analysis of |Skopal et al.|j2014[ >, who also showed 
that the spectral shape was a blackbody all the way out to the mid- 
IR. We confirmed our blackbody assumption for VI324 Sco by ex¬ 
tending our spectral fit up to the near-UV and comparing with flux 
values from Swift UVOT (Page et al. 2012^ . After applying redden¬ 
ing corrections from |Brown et al. ( 20 10^ our blackbody fit predicts 
a near-UV flux of 2.1 ± 0.4 x lO"*** erg s"' cm"^, a value in agree¬ 
ment with the observed flux, 2.4±0.16x 10"*^ erg s"' cm"^. As the 
near-UV flux is a factor ~ 10^ times lower than in the optical, this 
excellent agreement validates the assumed blackbody spectral pro¬ 
file. Results for the total optical/NIR flux and best-fit temperature 
as a function of time are shown in Figure]^ 
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As only three nights have complete NIR + BVR color data 
in nova Sco, we also consider separately the flux ratio on other 
nights with LAT detections obtained using just the V magnitude 
to calculate the optical/NIR flux, assuming a bolometric correction 
identical to that measured with better color data on day 15. Uncer¬ 
tainty in the total flux is estimated by combining the uncertainty in 
the best-fit constant olTset with an uncertainty in the temperature fit 
derived using the error in individual flux measurement. 


3.2 Gamma-Ray 


The Fermi LAT data for nova Sco and Del are taken from I Acker- 1 
|mann et^j2014] l. Spectral fits to the gamma-ray emission (over 
its entire duration) are provided by these authors for two profiles: 
Power Law (PL, dNyjdE oc and Power Law with Exponen¬ 
tial Cut-Off at energy E = (EPL, dNyjdE oc £-•' exp-^/^0- The 
best-fit values of the free parameters, T (PL) or s and Ec (EPL), are 
provided along with their uncertainties. 

Also provided for each day is the average photon number 
flux above 100 MeV, Ny oc ‘^dE. The energy flux Ey oc 

iiooMeV calculated from Fy according to 


Fy = 


N 


r EEhyi 

r JlOOMeV ^ dE “ 
^100 GeV dNy 
JlOOMeV 


(24) 


where uncertainties in Fy are derived from both the quoted uncer¬ 
tainties in Ny and in the best-fit parameters of the spectral fits. Al¬ 
though we calculate the gamma-ray flux each day, our analysis as¬ 
sumes a constant spectral shape over the duration of the outburst, 
as necessitated by the poor statistics of the gamma-ray detections. 


3.3 Flux Ratio 


Our results for the flux ratio L^/Lopt = fy/Fopt as a function of 
time for nova Sco and Del are shown in Eigure Key results are 
summarized in Table The ratio Ly/Lopt is consistent with being 
constant in time in nova Sco, but nova Del shows evidence for a 
moderate secular increase. 

We find an average value Ly/Lopt ~ 10“^'^(10“^'^) for nova 
Sco(Del), respectively. Assuming radiative shocks and gamma-ray 


efficiency Cy < 0.2 (5 2.3.1 2.3.2 1 , the minimum non-thermal par¬ 
ticle acceleration efficiency corresponds to > 0.1 - 0.01 in 
nova Sco and fnth > 10“^ in nova Del. This greatly exceeds the 
acceleration efficiency expecfed for electron acceleration at non- 
relativistic shocks, disfavoring lepton scenarios for the gamma-ray 
emission. Eurthermore, assuming a proton acceleration efficiency 
of Lith ~ 0.01 - 0.1, we conclude fhat a fraction /si, > 10 - 100 per 
cent of the optical emission was shock-powered in nova Sco. The 
shock-powered fraction is /h > 1-10 per cent in nova Del. 


4 DISCUSSION AND CONCLUSIONS 

The high densities within the nova ejecta at early times coincident 
with the observed gamma-rays imply that the shocks responsible 
for this emission are likely to be radiative ( §2.1| (. Though not ob¬ 
served directly, thermal X-rays from the shocks are a bsorbed by 
the ejecta and reprocessed to optical frequencies (^2.2|p]Based on 


Table 1. y-ray/optical flux ratio in classical novae and implications. 


Nova 

< log 

^7 

^-opl 

> 'Pg^ntLin 

Ab) 

hh 

VI324 Sco 

-2.2±0.4 

-1.5±0.4 

0.13(e„,h/0.1)-> 

V339 Del 

-3.5±0.2 

-2.8±0.2 

0.016(e„tt,/0.1)-‘ 


Minimum fraction of shock powered placed into non-thermal particles (eq. 
calculated assuming radiative shocks and (conservatively) €y = 0.2 (^2.3.1| 2.3.2|. 

Fraction of optical luminosity powered by shocks (eq. |l3| ), 
calculated assuming radiative shocks and (conservatively) €y = 0.2. 




Figure 2. Total optical flux {top panel) and best-fit blackbody temperature 
{bottom panel) as a function of time since outburst for novae VI324 Sco 
{blue) and V339 Del {red). Time is measured starting on June 1, 2012 for 
Nova Sco (beginning of optical outburst), and starting on August 16, 2013 
for Nova Del (epoch of first gamma-ray detection, within days of the optical 
rise). 


the observed values of L-^/Lopt, for the highest gamma-ray efficien¬ 
cies 6nth ~ 0.1 expected in hadronic scenarios, we conclude that 
a significant fraction /sh > 10 per cent and > 1 per cent of nova 
optical emission in nova Sco and nova Del, respectively, is pow- 


^ Nova Del did show X-ray emission at later times, after the epoch of peak had expanded or became sufficiently photo-ionized to reduce the column of 

gamma-ray and optical emission. This presumably occurred after the ejecta neutral absorbing gas (see below). 
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t (days since outburst) 

Figure 3. Ratio of gamma-ray and optical luminosities, Ly/Lapt, for novae 
VI324 Sco (top panel) and V339 Del (bottom panel) as a function of time 
since outburst. In Nova Sco, large squares show nights with full NIR+VBR 
coverage, while small triangles show the flux ratio estimated using just the 
V-band magnitude adopting the same bolometric correction as measured 
with the full data color data on Day 15. 

ered by shocks. These limits approach unity for more realistic ef¬ 
ficiencies e^th < 0.01 expected from the difficulty of accelerating 
protons given the expected geometry of the magnetic field in the 
pre-shocked ejecta (see below). 

This shock-powered emission mechanism is similar to that 
at work in some core collapse supernovae when the ejecta 
from the exploding star interacts with dense circumstellar matter 
(e.g. |Chevalier & Fransson|1994[|Larsson et al.|2011^ , which in ex¬ 
treme cases powers the most luminous supernovae yet discovered 
( [Smith et al.|2007l ). Our finding that a large fraction of nova opti¬ 
cal emission is powered by shocks requires a revision of standard 
models that instead assume nova thermal emission results from the 
direct outwards transport of thermal energy released in the white 
dwarf envelope by nuclear burning (e.g. [Hillman et al.|2014| ). The 
ultimate energy source (the thermonuclear runaway) is not in ques¬ 
tion, hut shock-powered emission taps into the differential kinetic 
energy of the outflow instead of its thermal content advected or 
diffused from small radii. As suggested by [Metzger et al.[ ( [2014[ l, 
shock-powered optical emission might also help explain some of 
the irregular light curve shapes, including plateaus and secondary 
maxima, observed in some novae (e.g.,[Strope et al.|2010^. A sim¬ 


ilar mechanism of internal shocks within a time-variable outflow 
may power an optical transient from the remnant produced hy the 
merger of binary white dwarfs jBeloborodov[201^ . 

If the majority of nova optical light is powered by radiative 
shocks, then the ratio of optical and gamma-ray fluxes should he 
relatively constant in time, assuming that the microphysical pa¬ 
rameters of the shocks remain constant (and that the conditions 
Tp_p > V 4 /C and tcooi/f < 1 remain satisfied in hadronic and lep- 
tonic scenarios, respectively). Nova Sco is consistent with a tem¬ 
porally constant ratio Ly/Lopt, while nova Del shows evidence for 
the gamma-ray emission becoming relatively stronger with time 
(Fig. 0 - The apparent increase of Ly/L„pi in nova Del could result 
from the reverse shock thermal emission becoming less radiatively 
efficient with time, i.e. x ^ Q)- Alternatively, it may reflect 

the early transition in this event to a nebular phase (as occurred on 
Day ll; [Skopal et al.|2014[ l, resulting in a loss of flux out of the 
measured optical band and into the UV. If nova shocks are indeed 
radiative, then we predict that Ly/L^pt should not exceed ~ 0.06 at 
any time, since this represents the limit of shock-dominated emis¬ 
sion /sh = 1 (eq. jl3| ) for realistic maximum efficiencies Sy < 0.3, 
e„,h < 0.2. 

The gamma-ray fo optical flux ratio also places a lower limit 
on the acceleration efficiency 6 , 41 , of relativistic non-thermal par¬ 
ticles at non-relativistic shocks. For nova Sco and nova Del we 
constrain e^th > 10“^ and 10“^, respectively (Table |^. Unlike tra¬ 
ditional studies modeling particle acceleration in supernova rem¬ 
nants, nova shocks evolve in real-time, in principle allowing for 
the study of time-dependent effects. Models used to constrain par¬ 
ticle acceleration in supernova remnants also sometimes require as¬ 
sumptions about the escape fraction of relativistic particles, e.g. if 
gamma-rays are produced by the collision of relativistic protons 
with nearby molecular clouds of irregular geometry. By contrast, 
nova ejecta serves as a relatively efficient hadronic “calorimeter” 
since relativistic protons are probably trapped when gamma-rays 
are observed (eq. GB)- 

Neither hadronic nor leptonic scenarios for nova gamma-ray 
emission can be ruled out by modeling of the y-ray spectrum alone 
( [Ackermann et al.|2014t . However, the tension between the value 
filth > 10 ^ we find is required in nova Sco and the much lower elec¬ 
tron acceleration efficiency e„th < 10 “^ inferred from observations 
of supernova remnants (e. g.[Morlino & Caprioli[20T^ and theoret¬ 
ical modeling jKato|20r4l[Park et al.[2014[ l appears to disfavor the 
leptonic scenario. Leptonic models also require a flat injected spec¬ 
trum p = 2 to be energetically feasible (eq. |23|). In contrast, ion 
acceleration efficiencies as high as ~ 0.1 are found for non- 
relativistic shocks by hybrid PIC simulations, but only in cases 
when the upstream magnetic field is parallel to the shock normal 
( [Caprioli & Spitkovsky[2014^ ; for perpendicular fields essentially 
no ion acceleration is seen, also inconsistent with our lower lim¬ 
its on e,ith' Hadronic scenarios involving ion acceleration at quasi¬ 
parallel shocks thus appear to be the most viable source of relativis¬ 
tic particles responsible for nova gamma-rays. 

However, this conclusion is problematic. Most of the power 
dissipated in the system occurs at the reverse shock, making it the 
most natural location for particle acceleration. If the fast outflow is 
ejected over a timescale which is long compared to the dynamical 
time near the surface of the hydrostatic white dwarf envelope (i.e., a 
quasi-steady wind), then its magnetic field will become dominated 
by its azimuthal component at the much larger radii where the re¬ 
verse shock occurs ( §2.3.l[ l. Such a field, being perpendicular to 
the radial shock normal, is not conducive to ion acceleration. 

This mystery could be resolved if particle acceleration oc- 
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curs at localized regions of quasi-parallel shocks, creating spe¬ 
cial regions of efficient acceleration ~ 0.1) with enth = 0 
across the larger bulk, thereby resulting in an average value of 
fnth ~ 10“^ - 10“^ consistent with observations. These localized 
regions of conducive field geometry could be caused by global 
asymmetries in the ejecta, such as oblique shocks occurring where 
the edge of the slow ejecta torus meets the faster bipolar wind 
(Fig. Alternatively, inhomogeneity of the nova ejecta (“dumpi¬ 
ness”) could also result in localized regions of oblique shocks be¬ 
tween the clumps. Indeed, low ejecta filling factors are inferred 
observationally by emission line modeling of the nebular phase 
(e.g., [Shore et al.|2013^ . Radial fields may also be produced near 
the shock as the result of Rayleigh-Taylor instabilities at the inter¬ 
face between the fast now outflow and the dense central shell (as 
may occur in supernova remnants; e.g., |Blondin & Ellison|200T| . 

Leptonic versus hadronic models for nova gamma-ray emis¬ 
sion may be further distinguished by their predictions for non- 
thermal X-ray emission. In nova Sco, Swift XRT observations 
placed an upper limit of Fx < 9 x 10“*'* erg s“* cm“^ on the 
0.3-10 keV flux on days 22-41 ( [Page et ^|2012^ , i.e. three or¬ 
ders of magnitude less than the simultaneous LAT-detected flux of 
Fy ~ 2 X 10“'° erg cm“^ s“*. Such a low value of FxIFy would 
at first appear to be barely compatible with the leptonic scenario, 
even when considering the most conservative case of p = 2 and 
a slow cooling spectrum vF^ oc from keV to ^ GeV energies 
(eq. (22|). The hadronic scenario fares even worse: secondary e* 
pairs from n* decay carry a similar total energy to that released in 
y-rays from decay. These pairs’ energies, ~ 0.1 - 1 GeV, corre¬ 
spond to Lorentz factors ~ 100 - 1000 that will upscatter ~ eV 
optical to energies ~ 10- 1000 keV with reasonably high efficiency 
(eq. |20|). However, a more detailed calculation must account for 
the important role of coulomb losses on the energy of the X-ray 
emitting electrons, which are likely severe due to the high densities 
of nova shocks; this will likely suppress the X-ray emission be¬ 
low that obtained by a naive extension of the IC or bremsstrahllung 
spectrum to lower energies. 

This dearth of early ~ keV emission from nova shocks could 
result from the high column of neutral. X-ray absorbing gas ahead 
of the shock at times coincident with the gamma-ray emission 
(Metzger et al.|2014 1 . However, higher energy X-rays > 10 - 100 
keV could escape without photoelectric attenuatioij^ motivating 
one to consider the prospects for nova detection with NuSTAR 
( [Harrison et al.|[2013| >, which has unprecedented X-ray sensitiv¬ 
ity above 10 keV. If a fraction fx of the LAT luminosity is ra¬ 
diated as X-rays of energy ex ~ 30 keV, this results in a num¬ 
ber flux of Nx ~ fxFyA,slex ~ 4 x 10“‘'(/x/0.01) s-‘ for a 
typical value Fy ~ 10“*° erg cm“^ s“* ( Ackermann_et__alj|20]_4^ 


and the NuSTAR effective area 200 cm'* jHarrison et al. 


[2013[ their Fig. 2). Given the NuSTAR estimated background of 
Vo- ~ 2 X 10“° counts s“* in this energy range across the point 
spread function jHarrison et al.[[^13[ their Fig. 11), we estimate 
that the integration time required for a 3cr detection is approxi¬ 
mately 9Na-INx ~ 100(/a:/0.01)“° ks. Reasonable constraints on 
fx ~ 0.01 could thus be achieved for a modest 100 ks exposure. 
The hard X-ray imager on ASTRO-H should also prove useful for 
nova follow-up. 

[Orio et ah] ( [2014^ reported NuSTAR observations of the fast 
symbiotic nova V745 Sco within 10 days of the outburst, placing 


° Although attenuation by inelastic electron scattering could be important 
for sufficiently high optical depths, leading to hard X-ray suppression. 


an upper limit of < 10“" ergs cm“° s“* on the > 50 keV emission. 
Gamma-ray emission was detected by LAT at 2 - 3cr significance 
and a flux of ~ 10“'° erg cm“° s“*, but lasting only for 2 days 
near the onset of the outburst ( [Cheung et al.|2014) , approximately 
a week before the NuSTAR observations and when (perhaps not co¬ 
incidentally in light of our results on shock-powered optical emis¬ 
sion) the optical light curve was at its maximum. By the time of the 
NuSTAR observations, the optical flux had decreased by a factor 
of ~100 from its peak value, such that if FyjFapt had remained ap¬ 
proximately constant in time (as in the novae in our sample; Fig.[^, 
then the NuSTAR upper limit correponds to fx < \0 in the notation 
above, i.e. not particularly constraining. We strongly encourage fu¬ 
ture NuSTAR observations of gamma-ray novae, ideally closer to 
the optical peak and coincident with LAT detections. 

In this paper we have focused on internal shocks in classical 
novae instead of symbiotic novae, because in the symbiotic case 
shocks and particle acceleration can result also from the interaction 
between the nova ejecta and the effectively stationary dense wind of 
the red giant companion jAbdo et al.[2010[|Martin & Dubus[2013[ (. 
Our analysis relies on the assumption of radiative shocks internal 
to the nova ejecta, as is justified by the higher mass outflow rates 
M ~ 10“°Mq yr“* and densities that characterize nova ejecta as 
compared to the lower values < 10“° Mq yr“* characteristic of red 
giant winds. Nevertheless, much of the same physics of relativistic 
particle acceleration could apply to both scenarios. Likewise, the 
presence of a dense equatorial structure associated with the nova 
ejecta could also be present in symbiotic systems (as was inferred 
in V407 Cyg; [Martin & Dubus|2013) . 

An additional prediction of the favored hadronic model is a 
GeV neutrino flux comparable to the gamma-ray flux that may be 
detected by future experiments ( [Razzaque et al.|2010^ . Depending 
on the maximum ion energy, novae may be important sources of 
TeV gamma-ray emission. [Sitarek et al.[ ( [20T5) recently presented 
upper limits on TeV emission from V407 Cyg by the MAGIC tele¬ 
scope ( [Aliu et al.|2012[ l which lie a factor of ~ 10 below the ex¬ 
trapolated Fermi LAT spectrum; similar TeV upper limits on V407 
Cyg were reported by VERITAS ( [Aliu et al.|2012^ . The observed 
cut-off in the photon spectrum may be due to an intrinsic cut-off in 
the accelerated ion spectrum (eqs. |16[ [18| ), or it may result from 
j-y absorption due to e- pair-creation attenuation by the nova op¬ 
tical light. TeV photons can pair create off fargef photons of energy 
£opt ^ 1 eV, resulting in an optical depth 


' nyf7"y_yRej 


: 0.3t. 


opt 


^opt 


^opt 


10°° ergs“ 


(25) 


near unity, where cTy-y ~ 10“°° cm° is the photon-photon absorp¬ 
tion cross section near threshold and riy = Lop(Topt/(4;rcR° £opt) is 
the energy density of the target optical photons near the shock (for 
Topt > !)■ We strongly encourage future additional TeV follow-up 
of novae near and after optical peak, even if not first detected by 
Fermi, as detections or upper limits can be used to place meaning¬ 
ful constraints on the maximum accelerated particle energy and the 
location of the gamma-ray emission. 

The sample of gamma-ray novae should expand rapidly in the 
next few years thanks to anticipated enhancements in the sensitivity 
of Fermi LAT resulting from improvements in its ability to perform 
low-level event reconstruction. The method for jointly analyzing 
optical and gamma-ray data employed in this paper will thus soon 
be applicable to a greater sample of events. Our results highlight 
the importance of obtaining broad (ideally bolometric) frequency 
coverage of nova light curves, including near infrared and ultravio¬ 
let wavelengths, at times coincident with the gamma-ray emission. 
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The utility of our model rests heavily on the assumption of radiative 
shocks, which depends sensitively on the quantity of mass and ve¬ 
locity of the nova ejecta, and how it is partitioned between the fast 
and slow components. Additional radio monitoring of gamma-ray 
novae is thus also key to better constraining the ejecta mass from 
these systems and its geometry, thereby strengthening or refuting 
the argument for radiative shocks. 
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